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We report on various nanocarbons formed from a unique structural pattern containing 
two pentagons, three hexagons and two heptagons, resulting from local rearrange- 
ments around a divacancy in pristine graphene or nanotubes. This defect can be 
inserted in sheets or tubes either individually or as extended defect lines. Sheets 
or tubes containing only this defect as a pattern can also be obtained. These fully 
defective sheets, and most of the tubes, present a very pronounced rippled (wavy) 
structure and their energies are lower than other structures based on pentagons and 
heptagons published so far. Another particularity of these rippled carbon sheets 
is their ability to fold themselves into a two-dimensional porous network of inter- 
connected tubes upon heat treatment as shown by hybrid Monte Carlo simulations. 
Finally, contrary to the common belief that pentagon/heptagon based structures are 
metallic, this work shows that this defect pattern should give rise to semi-metallic 
conduction. 
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I. INTRODUCTION 



Graphene may contain defects. This is now an established fact as it has been directly 
imaged, thanks to atomic resolution electron microscopy^} Defects in graphene have first 
been obtained by reconstruction around vacancies formed during electron irradiation and 
leading to pseudo amorphous domains made of pentagons, hexagons and heptagon d 1 * 2 * 4 !. 
They can also show some ordering as observed by Lahiri et alP, who describe an extended 
defect line, made of a periodic repetition of two adjacent pentagons and an octagon, in a 
graphene sheet epitaxially grown on a Ni substrate with a tilt boundary 

Defects affect the properties of graphene, for the bad in most cases but also sometimes for 
the good as in the case of Lahiri et al. where the line defect gives to the sheet the electronic 
properties of a one-dimensional metallic wir d 3 * 5 * 6 l The quest for stable space-filling one- 
or two-dimensional structural patterns has started long ago in the virtual world, actually 
much before the experimental characterization of such defects, thanks to ab-initio and force- 
fields based modeling approaches'^. Somehow, the computer-aided design of new graphene 
allotropes with targeted electronic or mechanical properties is even becoming a new area 
in carbon science^. Among all the structures proposed so far, hexagonal haeckelite H 5 e j 
has the lowest energy above graphene^. This can be explained by the absence of adjacent 
pentagonal rings and its high density of pentagon- heptagon (C 5 /C 7 ) pairs, a remarkably 
stable patterrPHEl. 

Two years ago, some of us have developed a new computational method based on a 
combination of image processing and atomistic simulation, the image guided atomistic re- 
construction (IGAR) method^, aiming at building atomistic models of dense, disordered 
yet nanostructured, carbons on the basis of their high resolution transmission electron mi- 
croscopy (HRTEM) lattice fringe images. In this approach, a 3D analogue of the 2D HRTEM 
image is constructed and serves as an external potential field, pushing the atoms towards the 
dark areas of the image during a liquid quench molecular dynamics simulation of a carbon 
system, in which carbon-carbon interactions are taken into account by a reactive empirical 
bond order potentia l 17 1 18 1 . 

A snapshot of a carbon sheet taken from a model of a rough laminar pyrolytic carbon^, 
reconstructed with this method, is shown Figure [T^i. As can be seen, the sheet is made 
of disoriented purely hexagonal domains (light grey) bound together by defective planar 
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FIG. 1. Snapshot of a carbon sheet taken from an atomistic model of a pyrolytic carbon, recon- 
structed using the IGAR method-^ (a). Bonds between threefold atoms belonging only to hexagons 
are shown in light grey, other bonds in dark grey; two (blue) lines show the dislocations, made 
of a succession of pentagon/heptagon pairs, between disoriented (as indicated by the (red) dot- 
ted arrows) hexagonal domains. The dashed red rectangle shows a space-filling structural pattern 
containing two pentagons, three hexagons (including a 90° rotated one) and two heptagons. This 
pattern can be used to form a fully defective carbon sheet (b), a line (c) or a point (d) defect in a 
graphene sheet. Dashed rectangles on panels (c) and (d) show the geometrical equivalence between 
the defect unit cell and the one of graphene. It is obtained from graphene by creating a divacancy 
(blue circles on panel (d)) and by a vertical to horizontal flip of the first neighbor parallel C-C 
bonds (green sticks on panel d). 

areas made essentially of C5/C7 pairs (dark grey) and more precisely lines of C5/C7 pairs as 
indicated by the blue solid lines. This agrees very well with some recent theoretical studies 
on dislocations in graphene showing the high stability of the C5/C7 patterrP^HUl Also, and 
even though the way this model was obtained has nothing in common with the irradiation 
of a single graphene layer, its structure is amazingly close from those reported by Meyer 



and co-workerd^. Looking at this figure, the structural pattern highlighted by the dashed 
rectangle appeared to be of a particular interest to us. Indeed, this domain, made of two 
pentagons, three hexagons (including a 90° rotated one) and two heptagons, allows for the 
construction of an entire carbon sheet based on this defect (see |l|b), noted ClClC^ hereafter, 
as well as for the insertion of a line (|T]c) or a point ([l]d) defect in graphene. Note that the 
latter point defect has been very recently observed in a HRTEMS1. 

In what follows, we report on the structure and properties of different nano-sheets and 
tubes based on this pattern as obtained from density functional theory (DFT) and force 
field calculations. Methodological details are given in section [TT| results are presented and 
discussed in section IIHI 

II. COMPUTATIONAL METHODS 

A. Density functional theory calculations 

Geometry optimization of the C\C\C* unit cell has been achieved using periodic den- 
sity functional theory calculations. The VASP package^^ was used with projector aug- 
mented wave pseudopotentials^^H to describe the valence-ion core interactions, and the 
Perdew-Wang-Ernzerhof generalized gradient approximation for the exchange-correlation 
potential. Atomic positions and lattice parameters were relaxed, and energy criteria of 
10~ 4 and 10 -3 eV/cell were used for electronic and structural convergences, respectively. A 
15 A interlayer distance was kept fixed for all calculations, and up to 450-point 2D grids 
were used for Brillouin zone sampling. The kinetic energy cutoff for plane-wave expansions 
was 400 eV. 

B. Hybrid Monte Carlo simulations 

Geometry optimizations of large C\C\C* planes and tubes were achieved using hybrid 
Monte Carlo (HMC) simulations^ coupled to a simulated annealing approach. The interac- 
tion energy between carbon atoms was described by the reactive AIREBO potential after 
removal of the van der Waals terms of this potential. Sheets and tubes were placed in an 
orthorhombic simulation cell with periodic boundary conditions in all directions. Sheets 
(respectively tubes) were placed on the xy plane (resp. along the z axis) of the box and a 
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large dimension (500 A) was given to the z dimension (resp. x and y dimensions) to ensure 
the simulation of isolated infinite systems. 

Two kinds of configuration updates were used (chosen at random) in our HMC implemen- 
tation, allowing for a proper sampling of configurational space in the isothermal-isobaric 
(NPT) ensemble: 

• Short constant energy molecular dynamics (MD) trajectories with initial velocities 
drawn from a Maxwell distribution at the suited temperature accepted with the usual 
canonical HMC acceptance criterion^. 

• Random contraction or dilation of the volume in one dimension of the system (x or 
y for sheets, z for tubes) accepted or rejected with the standard isothermal-isobaric 
Monte-Carlo acceptance criterion^ 

Simulations were initiated from the fully planar sheets or fully cylindrical tubes with cell 
dimensions taken from the corresponding graphene sheets or pristine tubes. Constant NPT 
HMC simulations of 10 4 steps were then performed at zero pressure and a given temperature 
(note that the position of one carbon atom with respect to the simulation cell was held fixed 
in order to avoid global translation or rotation of the system). After that, the temperature 
was scaled by a factor 0.9 and the process repeated until convergence of both potential 
energy and geometry. In the HMC simulations, the probabilities for MD and volume move 
attempts were respectively of 70 % and 30 %. The length of the MD trajectories (both in 
terms of the number of steps and of the timestep length) as well as the magnitude of the 
maximum volume change attempt were periodically adjusted to yield acceptance ratios of 
around 50 % for both kinds of moves. We used an initial temperature of 500 K and stopped 
the simulated annealing when the temperature was below 10~ 2 K, which appeared to be 
sufficient for a full convergence of geometry and energy. 

In the last part of the paper (see Figure 5) HMC was not used anymore to cool a system 
down to OK but to heat it up, as slowly as possible, to high temperatures. In that case 
similar simulations were used except that (i) the initial configuration was taken from the 
previously minimized system at K; that (ii) 5 x 10 4 HMC steps were performed at each 
temperature and that (iii) the target temperature was slowly increased by steps of 10 K. 
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III. RESULTS 



«— a=4.70lA -► 




FIG. 2. Optimized unit cell (a) and relative stability and potential well of the C^CqCj sheet as 
compared to graphene and hexagonal haeckelite (b) as computed from periodic DFT calculations. 
Red, blue and green labels on panel (a) respectively indicate unit cell parameters (in A), bond 
lengths (in A) and angles (in degrees); red circles, blue squares and green crosses on panel (b) give 
the energy as a function of the relative cell parameter for C^CqCj, haeckelite and graphene (taken 
as the reference for energy and cell dimension) sheets respectively. 

Figure |2^i shows the unit cell of the Cf Cf Cf as optimized by periodic density functional 
theory (DFT) calculations. No significant deviation from the rectangular lattice and no 
departure from planarity were found in our DFT calculations; however, this holds for the 
14-atoms unit cell of Figure [2^, and does not exclude non-planar configurations in supercells, 
i.e. with larger degrees of freedom, as we will see in what follows. The energy of the 
optimized structure (see Figure^) is only of 227 meV/atom above graphene, making it 
the lowest energy defect-containing carbon sheet reported so far. To our knowledge, it was 
formerly the case of hexagonal haeckelite, with an energy of 240 meV/atom above graphene 
(note that we compare to those structures as optimized using the same methodology that 
we applied to the C\C\C?j sheet). Also, Figure^ shows that the potential wells of the three 
structures are very similar. 

We present now the properties of large carbon sheets and nanotubes based on the C\C\C% 
pattern as obtained from cell and geometry optimizations using hybrid Monte CarlcP (HMC) 
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FIG. 3. Snapshots of a periodic 5 x 21 CfCgCf sheet (a,b), a 1 x 21 line defect (c,d) and a 
single point defect (e,f) inserted in graphene sheets. Geometries have been obtimized using HMC 
simulated annealing and the AIREBO force field. (a,c,e): plane view; (b,d,f): lateral view. 

simulated annealing and the AIREBO force flelcP^. Figure [3] shows a front (a) and a side 
(b) view of an optimized periodic CfCfC| sheet made of 5 x 21 unit cells. The energy of 
this sheet, as given by the AIREBO potential is of —7.131 eV/atom, 279 meV/atom above 
graphene as computed with the same potential (—7.408 eV/atom). Relative energies of 
the C^CqCj, haeckelite and graphene sheets given by this potential are in fair qualitative 
agreement with DFT calculations (haeckelite is found at 323 meV/atom above graphene 
with this potential). However, a major difference is observed between haeckelite and Cf Cf C| 
sheets. Indeed, if the geometry optimization of a large haeckelite sheet reveals an almost 
planar geometry, Figures [3^i and[3]D unambigously unravel a wavy structure developing itself 
along the a unit cell direction, the direction of the adjacent hexagons lines in the structure. 
Optimizing the number of unit cells in that direction we found out that the best possible 
periodicity is four unit cells along the a axis, giving rise to an optimal energy of —7.135 
eV/atom. The line defect of Figures^ and|3}i is also of particular interest. First, and even 
though it is not easy to directly compare the energies of different line defects, it should be 

7 



of rather low energy as it is made of a remarkably stable pattern (the CfCfCf). Second, 
while other observecP or predicted^ line defects usually induce a curvature of the sheet 
perpendicularly to the line, Figure [3]i clearly shows that the C\C\C^ pattern gives a wavy 
structure to the sheet, along the line defect, in close similarity with the wavy structure 
developed on the perfect Cf Cf Cf sheet, though attenuated by the large hexagonal domains 
on both sides. The point defect of Figures^ and[3]F also has an excellent energetic stability, 
but affects only very locally the structure of the sheet. It is an antisymmetric defect with 
an inversion center localized on the center of the inverted hexagon. We recall that this point 
defect has been very recently imaged by Kotakoski et al. during a HRTEM observation of 
an irradiated graphene layer^. 




□ □ 

FIG. 4. Snapshots of carbon nanotubes built from the Cf CqC| . Geometries have been obtimized 
using HMC simulated annealing and the AIREBO force field. In analogy with pristine nanotubes, 
tubes having the lines of hexagons along (b,c) and perpendicular to (a) the tube axis are respectively 
noted as zigzag (ZZ) and armchair (AC), a: side view of a zigzag tube made by folding a ribbon of 
nine unit cells along the a axis (ZZ9); b: side view of an armchair tube made by folding a ribbon 
of five unit cells along the b axis (AC5); c: side view of an armchair tube made by folding a ribbon 
of six unit cells along the b axis (AC6). (d,e,f): top views of tubes (a,b,c). 

As already said, the Cf Cf C| pattern likes curvature. This should then be a particurlarly 
stable pattern to build nanotubes with. In close similarity with pristine, graphene based 
zigzag (ZZ) and armchair (AC) tubes we have built CfC|Cf ZZ and AC tubes by folding 
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FIG. 5. Energy of zigzag tubes as a function of their diameter as obtained from hybrid Monte 
Carlo optimizations (the origin of energy is the lowest energy CfCfCy sheet). 




o 
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FIG. 6. Same as Figure [5] for AC tubes. 
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C^CqCj ribbons along respectively the a and b axes of the pattern's unit cell (see also the 
caption of Figure |1J). We present now some results of HMC geometry optimizations on ZZ 
and AC tubes of different diameters, ranging from 7 to 30 A, and with lengths around 9 nm. 
Figure [4] present snapshots of the optimized ZZ9, AC5 and AC6 tubes. Their diameters, 
taken as the double of the mean atomic distance from the tube axis, are respectively 13.2, 
12.7 and 14.7 A. All these tubes have rather low energies, respectively —7.141, —7.084 and 
—7.101 eV/atom for ZZ9, AC5 and AC6; and present different morphologies. All zigzag 
tubes studied in this work have a "necklace of pearls" shape as shown on Figures |4^i and 
|4jd for ZZ9, in close similarity with some haeckelite tubes^. Low diameter armchair tubes 
(from AC3 to AC5) have a fairly cylindrical morphology (see the snapshots of tube AC5 on 
Figs [4|3 and |4^), typical of graphene based tubes. However, from AC6 (14.7 A) to AC12 
(29.9 A), armchair tubes develop a pronouced helical morphology (see the snapshots of tube 
AC6 tube on Figs^ and[|f). 




FIG. 7. Evolution of the energy with respect to temperature during a HMC simulation of the 
heat treatment of a periodic 17 x 9 C|C|Cy sheet from to 2000 K (empty circles) and cooling it 
back to K (filled circle). The same simulation as applied to a haeckelite sheet of similar size is 
indicated shown for comparison (empty squares). 

On Figure [5] we plot the energy of ZZ tubes as a function of their diameter relatively to 
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the sheet of lowest energy (-7.135 eV/atom). As expected, the energy per atom, around 40 
meV above the C\C\C^ sheet for the smallest tube, diminishes with the tube diameter at 
low diameters thanks to the relaxation of curvature stresses. However, the behaviour of the 
curve above 10 A is particularly unusual. First, we can see that ZZ tubes with diameters 
in the range [12:25 A], from ZZ8 to ZZ13 here, have lower energies than the rippled sheet. 
Second, the curve goes through a minimum around 15 A, with the minimum energy tube 
found here, ZZ11, being 7 meV/atom lower than the sheet. Finally, above 25 A diameter, 
ZZ18 tube has a positive energy with respect to the sheet. These results indicate once again 
the high affinity of the C\C\C^j pattern for curvature along the lines of hexagons (a axis) 
and tend to show that there exists an optimal curvature for this pattern which corresponds 
roughly to a 15 A diameter ZZ tube. 

Figure [6] is equivalent to Figure [5] for armchair tubes. The first points to note are that 
AC tubes present higher energies than ZZ tubes at equivalent diameters and that all the AC 
tubes studied in this work have higher energies than the lowest rippled sheet. For instance, 
the most stable AC tube studied in this work, AC10 (0 = 30 A), is found 20 meV/atom 
above the sheet. Unlike ZZ tubes, the evolution of AC tubes energy with respect to diameter 
is monotonically decreasing; however, a step is observed between the low diameter cylindrical 
tubes (AC3 to AC5) and the larger helical tubes (AC6 to AC10). 

In order to investigate the thermal stability of the C\C\C% sheet we now present results 
of a HMC simulation in which a periodic sheet of 2142 atoms was slowly heated from 
to 2000 K. The evolution of the energy per atom as a function of temperature is shown 
on Figure [7] and some snapshots of the system during the heating process are given Figure 
[8| As can be seen on Figure [7j the energy evolves rather linearly with temperature at low 
temperatures and the sheet shows three well defined waves (see Figure^). Around 1000 
K, a slight drop in the energy increase is observed and the sheet evolves from three low 
amplitude waves (Figure [8^,) to two well-pronounced waves (Figure^). The energy then 
keeps on increasing linearly until another drop at around 1400 K. At that point, two parts of 
the sheet enter into contact and form a regular bridge, containing sp 3 hybridized atoms, all 
along the b unit cell axis (Figure ^p)- This is shortly followed by the formation of a second, 
and symmetric, bridge around 1700 K to form an almost perfect periodic two-dimensional 
network of interconnected tubes (Figure [8ji). This rippled nanocarbon obtained at 2000 
K, shown on Figure [8jd., has then been optimized by a HMC simulated annealing quench 
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FIG. 8. Snapshot of the Cf Cf C? sheet of Figure @ after annealing at 500 K (a), 1000 K (b), 1500 
K (c) and 2000 K (d) . A snapshot of the system cooled back to K (filled circle on Figure [7]) is 
also shown (e). (note that five periodic replicas of the system, along the horizontal direction, are 
shown on (e) to highlight the 2D nature of the material and its network of aligned (ID) porous 
channels). 

down to K, the latter structure being shown Figure [8^. Its energy, of —7.167 eV/atom, is 
32 meV/atom below the one of the initial sheet, indicating, once again, the affinity of the 
C^CqCj pattern for high curvatures. The diameters of the interconnected tubes forming this 
porous structure are actually very similar to those of the ZZ9 tube presented on Figures [pi 
and|i]i, and correspond to the optimum curvature for zigzag tubes (see Figure [5]). Looking 
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FIG. 9. Enlargement of Figure [8]f with two coloured regions highlighting the open windows to an 
upper tube (red) and a lower tube (blue). 

at Figure [8^, the structure of this carbon immediately reminds of some shock protection 
materials used in everydays life packaging and there is no major risk taken in guessing that 
its mechanical properties would be particularly unusual. Also, it should be regarded as a 
very interesting adsorbent due to its aligned network of high aspect ratio tubular pores. 

Figure [9] shows an enlargement of Figure [8ji, to unravel the porous nature and accessibility 
of this new form of carbon. As can be seen on this figure, those pores are extremely easy 
to access, thanks to numerous windows all along the structure, probably facilitating any 
adsorption/desorption process with respect to other CNT based adsorbents. 

Finally, it is interesting to note, as illustrated on the top of Figure]?], that the same heating 
simulation on haeckelite shows a perfectly linear evolution of energy with temperature which 
corresponds to a perfect conservation of the sheet structure (not shown), aside from thermal 
vibrations. 

IV. CONCLUSION 

We have proposed new rippled carbon planes and tubes based on a new structural pattern, 
the G\C\C^ pattern, recently observed in some atomistic models of pyrocarbons. These 
structures can be either entirely made of this pattern or inserted as line or point defects in 
graphene or graphene-based tubes. Both DFT and force field energy minimisations show that 
structures built from this pattern have low energies (actually, the lowest energies reported 
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FIG. 10. Electronic density of states (DOS) of the C|CgC| (straight line) and graphene (dashed 
line) unit cells as computed from DFT calculations. 



so far) for structures with a high density of non-hexagonal rings. Another evidence of the 
high stability associated to this pattern is that the point defect of Figures [3^ and [3]F has 
been recently imaged in a TEIvP. Structures based on the C\C\C^ pattern love curvature 
and show rather unusual properties: (i) sheets have a pronounced wavy structure; (ii) line 
defects develop a parallel wavy structure; (hi) tubes may show lower energies than sheets; 
and (iv), sheets are able to close themselves under heat treatment to form a 2D ordered 
nanoporous carbon. Finally, the highest surprise stemming from this defect type maybe 
lies in its electronic properties. Indeed, it is now a common thought that pentagon and 
heptagon containing carbon structures are can be seen on the density 



of states plot of Figure [TOj and despite its high density of pentagons and heptagons, the 
C^CqCj network is semi-metallic, alike graphene (although its DOS has more features than 
the one of graphene due to its lower symmetry and the more complex connectivity path of 
its carbon network). 
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